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ABSTRACT 



We present the convolution integral for fermion pair production in the electroweak standard 
theory to order O(a) including also soft photon exponentiation. The result is complete in the sense 
that it includes initial and final state radiation and their interference. Prom the basic result - 
analytic formulae for the differential cross section - we also derive the corresponding expressions for 
total cross section errand integrated forward-backward asymmetry Aps- The numerical importance 
of different contributions for the analysis of experiments at LEP/SLC energies is discussed. 
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INTRODUCTION 



Stimulated by the possibility of very precise experiments at e + e - storage rings, the QED 
corrections for the reaction 

e+(k 2 ) + e-(h) — (7, Z) — f+( P2 ) + r( P i) + n 7 (p) (1) 



have been studied in very detail for a large energy region including the Z resonance. A comprehensive 
collection and comparison of available results within the electroweak standard theory may be 
found in o and in references quoted therein. Analytic and semi-analytic formulae proved to be of 
great value both for a deeper understanding of process ( [l]) and for ensuring a high reliability of 
Monte Carlo codes ® . Recently theanalytic approach to QED has also been used for the analysis 
of Z line shape data at LEP/SLC Q. The line shape is the integrated total cross section ap (s) 
as function of s = 4E 2 where E is the beam energy. This most inclusive observable has been 
studied theoretically since many years and is well-described to order 0(a ) B ^llOt with soft photon 
exponentiation and also for several massive intermediate bosons , higher order leading 

logarithmic corrections $±i< tL§ and the complete QED initial state radiation corrections to order 
0(a 2 ) A compilation of many important contributions is contained in . 

Only recently, similar analytic results have been published for the integrated forward-backward 
asymmetry Aps ■ The complete 0(a) convolution integral with soft photon exponentiation for 
Aps was obtained in S3 (see also 11 ) and the leading logarithmic approximation for initial state 
radiation including 0(a 2 ) terms in EH . A survey on results about Afb is contained in EH . 

Compared to ap and A p b , analytic results for the differential cross section dcr/dcosO are scarce. 
Earlier attempts are HH . Other distributions are treated in —\ . In pure QED, the first compact 
analytic expressions for hard bremsstrahlung corrections to the differential cross section have been 



derived in 1^3 . A formalism for leading logarithmic approximations may be found in fc^l , though 
without application to du/dcos©. Recently, semi- analytic formulae allowing for quite realistic 
cuts have been applied in ^H) . This article contains the systematic presentation of an analytic 
calculation of QED corrections to the angular distribution for reaction ( |l|) . Though computer 
codes relying on the present results have already been released and used for a ppl ications in LEP 
physics, very few of the material has been presented so far ; see also H HHT^l . In chapter 2 
we introduce the notation and describe in short two derivations of the analytic expressions for the 
angular distribution. Chapter 3 to 5 contain the basic result which consists of compact, explicit 
expressions for the soft and hard photon radiator functions in the convolution integral due to initial 
(chapter 3) and final (chapter 5) state radiation and their interference (chapter 4) to order 0(a) . 
They also include higher order corrections and some formulae for ap and Aps ■ Chapter 6 contains 
numerical results and conclusions. 



2. NOTATIONS AND PHASE SPACE PARAMETRISATION 



The explicit analytic formulae for the differential cross section proved to be too cumbersome 
to be described in detail here. For that reason, we choose a semi- analytic presentation with the 
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following notation: 



da 
dc 



E 

m,n=0,l 



E 

=T,FB 



^ Re[a A (s, s;m,n) R A (c;m,n)]. 



The scattering angle c = cos is defined between the produced fermion / + (/ = fj,, v, q) and the 
positron beam. The functions a° A {s, s; m, n) are reduced Born cross sections defined below and the 
corrections R A (c; m, n) can be expressed by semi-analytic formulae for soft and hard parts of QED 
bremsstrahlung radiator functions S A (c,e;m,n) and H A (v,c): 



R\(c; m, n) 



(s, s; m, n) 



e)H%(v,c)}. 



(3) 



In ( H|), we allow for a possible cut-off A ( 0< A < l-4m^/s) on the energy v of the emitted photon 
in units of the beam energy. The effective energy s' of the created fermion pair is s' = (1 — v)s. The 
sums over m,n in ( |2|) include photon (m = 0) and Z boson (m = 1) exchange. A generalisation 
to the case of additional vector bosons Z n ,n >1, is straightforward. Indices a = e,f,i are used 
for initial and final state radiation and their interference. The angular dependence of the Born 
cross section has been formally included into the initial state contribution a = e. The P-even and 
P-odd cross section parts carry index A = T and A = FB. Under CP-invariance, they are also 
C-even (C-odd), correspondingly . After symmetric (anti-symmetric) integration over cosG the 
functions {R%b' functions) yield the total cross section ctt (the integrated forward-backward 
asymmetry Aps )'■ 

da 



ax 



A 



Within our terminology, 



FB 



E 



ax 



1 

ax 



dc- , 
i dc 

1 , da 

dc 

o dc 



Re[dA a\ (s, s; m, n) R\(m, n)], 



» , da 
dc— 
i dc _ 



(4) 



(5) 



(6) 



m,n=0,l a=e,i,f 



R a A (m,n) 



d A l 



dc R\{c; m, n) 



A = T, FB, 



d F B = 1- 



(7) 



(8) 



The additional numerical factors in ( ||, [?]) ensure the usual total cross section normalisations 
for <tt and apB- 

The dependence of dc /dc on weak neutral couplings and electric charges of the fermions and 
on possible polarisations and weak loop corrections together with the typical Born-like resonance 



behaviour have been collected in the explicit Born factors a A ,0 (s, s; m, n). As a consequence, the 



Mi and width Tz = Ti 



13 and 



functions R A (c;m,n) depend only on the Z boson mass Mz 
on kinematic variables s, cos0, and the energy cut-off A. Of course, fermion masses which are 
assumed here to be small compared to s, Mz, Tz show up in certain logarithmic mass singularities. 

i A \ . /// . n, j have different energy dependence for a = e,i, f: 



The reduced Born cross sections ct°' (s, s'; m, n) 



a A (s, s ;m,n) = a A {s ,s;m, n), 
a A (s, s ;m,n) = a^s, s ; m, n), 



(9) 
(10) 
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a A ,0 (s, s'; m, n) = <r A (s, s; m, n). (11) 
Here, A = FB or T if A = T or FB, and 

°a (s,s';m,n) = ^-C A {m,n; Ai, A 2 , hi, h 2 )cf-\x m {s')x* n (s) + Xm(s)Xn( s ')]; ( 12 ) 

*„(*) = — L T' ( 13 ) 

= M2-iM„r n (s), (14) 
r «(*) * ^2 r - ( 15 ) 

lvl n 

In ( [l~2] ), the fermions have a color factor c/ = 3 in case of quarks and cj = 1 for leptons. Their 
vector and axial vector couplings Vf(n) and a/(n) to gauge boson Z n are contained in Ca- 

C T (m,n;X 1 ,X 2 ,h 1 ,h 2 ) = 

{Xi[v e {m)v* e {n) + a e (m)al(n)]+ X 2 [v e (m)a* e (n) + v* e (n)a e {m)\\ 
{hi[vf(m)v*f{n) + af{m)a*j(n)}+ h 2 [vf(m)a*f(n) + v j(n)a/(m)]}, (16) 

C FB (m,n;Xi,X 2 ,hi,h 2 ) = C T (m, n; A 2 , Ai, /i 2 , hi). (17) 

We allow for longitudinal polarisations A_, A + of both the electron and positron beams and for final 
states with helicities Due to the CP-invariance of the problem and the (v,a) structure of 

the interactions, the following combinations of polarisations are possible: 

Ai = l-A+A_, A 2 = A+-A_, (18) 

hi = ±(l-h+h-), h 2 = ~(h+-h-). (19) 
In case of the standard theory HI , we use the conventional couplings of photon and Z boson: 

go = e, v f (0) = Q f , a f (0) = 0, (20) 

51 = (v^Mf) 1 / 2 , v f (l) = I^(f)-2Q f sm 2 Q w , a f (l) = J 3 L (/), (21) 

where I^(f) is the third component of the weak isospin of fermion /; 2 I^(e) = Q e = -1. 

The inclusion of non-QED weak loop effects (see @ and references quoted therein) is 

trivial if one uses the form factor approach for their description . Coupling constants go,gi, and 
the vector couplings v e (l),Vf(l) become complex and s-dependent. Further details are explained in 



19] 



For an exact treatment, an additional form factor v e f(l) replaces the product v e (l)vf(l) in 
( |l(]) whenever it appears there. In principle, energy-dependent form factors should be understood 
to be part of the integrand in ( ||, ||) . We have checked by explicit comparisons that, due to 
the minor s-dependence of them, the numerical error which is implied by the approximation is 
practically negligible 0' 13 . 

For notational convenience, all angular dependences have been included into the radiator func- 
tions. The functions R\(c\ m, n) of ( |3|) are result of an incoherent sum of real and virtual photonic 
corrections. They are obtained by straightforward but lengthy Feynman diagram calculations with 
extensive use of the analytic manipulation programs SCHOONSCHIP @ and REDUCE 63. The 
Rji(c;m,n) are gauge-invariant and, if integrated without cuts, also Lorentz- invariant. Determi- 



nations of the vertex and box corrections "I , soft photon and hard photon contributions deserve 
different techniques, correspondingly. 
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We would like to give short comments on the two methods which have been used for the 
0(a) bremsstrahlung phase space integrations. In one approach , we used the so called R 7 - 
system, the rest system of (/ — 7), where the three-momentum relation p(f~) + pX'y) = is fulfilled. 
The phase space parametrisation chosen is: 



7T / c - m / s 

dT = — t / dcos© / xdx — 
As 6 J- Cm Jo s 



x + mj Air J -i 



1 fl r 2-K 

— \ dcosG^ 



(22) 



where (O^ , <p^) are photon angles in the i? 7 - system and x the momentum fraction of / + in the 
cms in units of the beam energy: 



x = -2p(f + )\p(e + ) +p(e~)] = 2^E(f 



(23) 



Some details of the calculation have been described in HH , includi ng t he n-dimensional treatment of 
soft photon contributions which follows the methods develope d in 153 . The hard photon integration 
has been done with SCHOONSCHIP using tables of integrals $3 . The total cross section orand the 



integrated forward-backward asymmetry without cut ff have been determined with this technique 
as well as the angular distribution which remained unpublished so far with exclusion of the pure 



QED case 



' . Another set of integration variables allows for a cut on the photon energy: 

7-2 rC m rVm rVlmam 1 rl-K 



As 



d cos Q dv 

J0 J D 2 , 



1 

-V 

2 



1 ± V m (S 



Am 2 f/s. 



,max(min) 

v m (s) = 1 

In ( [24|) , (p~ is the azimuthal angle of the photon in the cms and 

2 



dv 2 -J o d^, 



Ml/2 



V 2 



-p(/ + M7) = l-2£(/-)A/i. 



(24) 

(25) 
(26) 

(27) 



We should also comment on the integration boundaries for cos©: 

c m = 1 - 2ml/ s. 



(28) 



The radiators for the angular distribution contain terms behaving in the hard photon corner like 
ln(l±c), (l±c) -2 etc. which become singular and for initial state radiation even non-integrable at 
c = ±1. In fact these quantities are related to t = —2p(e + )p(f + ) and x: 



t 



x 
2 



1 



Amis 



2 — AmjrS cos 0, 



-[1 

2 L 



(1 - 2m 2 e / s) cos ©]• 



(29) 



(30) 



From ( |28|), it becomes evident that the phase space integrals ( 2^,24]) are properly regulated. In 
order to get a convolution representation for dcr/dcos B one has to perform in ( 24) two integrations. 
This has been done with SCHOONSCHIP and REDUCE, based on tables of integrals for hard 
bremsstrahlung. The soft photon part is the same as determined earlier; see, e.g., in 01131. While 
all the hard photon parts of the radiators for the angular distribution remain independent of the 
type of exchanged gauge boson, this is not the case for the box diagram contributions to the 
soft photon part of the QED interference corrections as has been indicated in the notation in 



1. 



( |3|) . After summing up the vertex or box contributions with those due to soft photons, the 
resulting S A (c, e; m, n) are infra-red finite but yet dependent on the infinitesimal parameter e which 
discriminates between soft and hard photons. The integral over the photon energy in ( |) as a 
whole is independent of e as may be seen from the formulae of the following section. Restricting 
oneself to a semi-analytic (i.e. containing one numerical integration) result for dcr/dc, one can take 
into account even cuts on the acollinearity and on the fermion energies as has been demonstrated 
recently in HH where also ( ^4]) was the starting point. 



3. INITIAL STATE RADIATION 



In this chapter, we give a systematic presentation of initial state corrections. 



3.1. THE RADIATOR FUNCTIONS FOR THE ANGULAR DISTRIBUTION 



There are two radiator functions for the QED inital state corrections, both of them containing 
well-known soft photon parts and hard photon parts which will deserve some comments. The soft 
photon parts including vertex corrections are: 



S A (c,e;m,n)=D A {c)[l + S(e,(3 e )}, 

D T (c) = 1 + c 2 , 
D FB (c) = 2c, 

where we include in ( also the Born cross section. Further, 



T,FB, 



(31) 

(32) 
(33) 



S(e,p e )=p e lne + - + -Q 



a 



7T 



7T 

T 



(34) 



e = — Q 2 e (L e -l) 

IT 



In -At. 

mi 



(35) 



The hard photon parts are symmetrised for the C-even part (A = T) and anti-symmetrised for the 
C-odd part (A = FB): 



Htfb(v,c) 



a 



7T 



Qe[hT.FB(v, C) ± h e TFB (v, -c)]. 



As already stated, the hard radiator parts depend on only two variables v, c: 



h F (v,c) = x 



L. 



2z 



2z 



10 



~^r 2 r 2 r 1 H + (-— r 4 + —r 2 - 4z) + - — (r 2 r 3 



7 



V 



1 



3z 3 
44 

^y3 



26 

^2 (3r 4 + 8zr 2 + —z 2 ) + -3 (8r 3 + ^-zn) - ^ ( ^fr 2 + 4z 



3jz 

2 '22 
3 



h e FB (v,c) = 



L c . 2z. 2 4 . , 2z . 
r2 -o (n ) + -»"2 ^^1(^2 + z) + ^-(3r 2 + 22; J 



(36) 



(37) 
(38) 
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The following abbreviations are used: 



7 2 

L c = ln^- + L e , (39) 

z=l-v, r n = l + z n , (40) 

c± = -(l±c), 7 = c+ + zc_. (41) 

As one should expect, the distributions are singular for v approaching 1 (hard photon emission) 
and c approaching ±1 (collinear radiation). In the soft photon corner (v — ► e), the hard radiators 
become singular too: 

]imH%(v,c) = ^D a (c) + 0(l), A = T, FB. (42) 

From ( |42| ) the necessary compensation of the dependence on the soft photon parameter e in ( ^J) is 
evident. Performing the integration over the photon energy variable v leads to analytic expressions 
for the initial state corrections R A (c;m,n) to du/dc as defined in ( §-||). These integrations are 
cumbersome but not too complicated. The structure of the integrands cr A °(s, s'; m, n)H A {v, c) 
leads to integrals containing rational functions, logarithms and Euler dilogarithms. The interested 
reader may envisage the Fortran program MUCUTCOS in the package ZBIZON E§l where the 
corresponding analytic expressions for the angular distribution are coded. We only would like to 
comment on the origin of the radiative tail within the present formalism. For initial state radiation, 



it is useful to linearise in s' the resonating function ( [12-13) 



777^ S 1 — 777.2* m 2* _ m 2 



1 1 



777^* s' — ml 



if at least one of the interfering gauge bosons is massive. For n = p, there arises the factor 

s i M n s 



ml* - ml 2 T n Ml ' 



(43) 



(44) 



setting the scale of tail effects, M n /T n . The threshold is defined by the onset of influence of the 
imaginary quantity ( ^4|) onto the real type cross section. So, it needs another imaginary quantity 
which is found from (s'-m^) -1 . After integration over the photon energy, functions of the following 
type arise: 

^ dv- l — = - In ]~* n ~l = ~L Rn (A), (45) 

1-v- R n 1 - R n - A 

Rn = m 2 Js. (46) 

If the cut-off A is not as infinitesimal as e (i.e. if hard photons are radiated), a substantial imaginary 
part is developed if s > M% . 

The photonic case, n = p = 0, must be treated with care. In order to get the exact hard 
bremsstrahlung correction for pure photon exchange one must use a reduced Born cross section 
( [l2]) which includes an additional phase space factor under the integral: 

a° A (s', s'; 0, 0) yjl - \m)js' (1 + 2m}/s') a° A (s', s'; 0, 0). (47) 

The additional threshold factor in ( |47]) influences only a minor part of the calculation and it is not 
too difficult (though lengthy) to take it into account in an analytic calculation. The pure photonic 
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corrections to the angular distribution may be found in eiS . The function Fq defined there in 
eqs.(5), (10) corresponds to R^(c; 0,0) introduced here in ( |3|). Interesting enough, in the existing 
literature there is no explicit mentioning of the correction R e FB {c; 0,0) or of its integral which 
contributes to Afb (see below). Since QED corrections due to photon exchange do not contain 
axial- type couplings, there does not arise any contribution of Rp B (c;0, 0) to the cross section - if 
not at least one beam polarisation and one final particle helicity are nonzero; see ( pit- 19). 



3.2. SOFT PHOTON EXPONENTIATION 



As may be seen from ( ^2j), the hard parts of the radiator functions become singular in the soft 
photon corner of the phase space. Their divergent part can be combined with the soft photon radi- 
ator ( 31,34) in order to get the lowest order of an infinite sum over soft photon contributions. The 
treatment can follow exactly the arguments given in 113 for the total cross section. The result 
is the following modification of the initial state radiative corrections for the angular distribution: 

R A (c;m,n) = f A dv a \ { f> ^ w) {D A (c)[l + S{(3 e )](3 e v^- 1 + H%(v, c)}, (48) 

H A (v,c)=H A (v,c)-^D A (c), (49) 

v 



S(J3 e ) = -Q 

7T 



^ 1 3 /r 



(50) 



In fig. 1 the contributions of the initial state radiation corrections to the differential cross section 
are shown as functions of the scattering angle. 

We choose as typical energies yfs = 30 GeV (TRISTAN); 91.1 GeV (the Z boson mass value; 
LEP,SLC) and 200 GeV (tail region, LEP 200). The relative importance of the cross section 
contributions depend on three different components; the coupling combinations C^m, n) which 
are in addition channel dependent, the factors |x| 2 ) which have been split away for reasons 
explained above, and then the Born plus QED correction factors. Shown are the contributions for 
soft photon exponentiation. At lower energies, pure photonic corrections naturally dominate for 
ot while they are zero for (Jfb (if there are no polarisation phenomena). At the Z peak and beyond 
in the tail region, the Z exchange dominates. For | cos 0| approaching 1 (beam directions), the cross 
section rises extremely due to the collinear hard photon emission. In fig. 2 the net Born plus initial 
state radiation cross section is shown as function of the scattering angle with the photon energy 
cut-off A as parameter. In the lower energy range the angular distribution is nearly symmetric due 
to the dominance of pure QED. The same is true at resonance, but here it is due to the smallness of 
the coupling combination accompanying the nonsymmetric pure Z exchange contribution (in case of 
muon production). In the tail region, a pronounced non-symmetric angular behaviour occurs. Hard 
bremsstrahlung leads to a rising of the cross section with exclusion of the resonance region where 
due to initial state radiation the reduced effective energy falls below the peaking value. A cut on 
the maximal photon energy reduces the cross section while soft photon exponentiation compensates 
at least partly the negative virtual plus soft photon terms. 



3.3 INTEGRATED CROSS SECTION AND FORWARD BACKWARD ASYMMETRY 
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Integrating over the angular dependence of the initial state radiation ( EHJ 36) , one gets the 
radiator functions for ot ^ and Afb '— • These are definitely different from each other already in 
the 0(a) leading logarithmic approximation. To some extent, this remained unclear in the comment 
on that point in lia. Nevertheless, at the Z peak the differences are small e3. For more details 
we refer to 113 and references therein. 

The one remaining integral over the photon momentum can also be calculated. For that purpose, 
one has to take care of the energy-dependence of the width function in the Z boson propagator ( |i~3| ). 
Sufficiently far away from fermion thresholds, approximation ( |l^) is excellent. This fact allows to 
apply the Z boson transformation H thus ensuring that the s-dependence of the Z width does not 
complicate the analytic integrations: 

X'i(s) = xi(s), (51) 

IV ^ - 1 



xito 



47TO 



1 + i 



M z 



M' Z = M Z U1 + 



Ml 



M'z 2 + iM' z T' z ' 



M z + A z , 



A; 



1 7 



2M Z 



-37MeV, 



(52) 

(53) 
(54) 
(55) 



where we used M z = 91.1 GeV and T z = 2.6 GeV @ . As a result, one gets for the total cross 
section ( |j) the QED correction functions E^(m,n) to order 0(a) and E^(m,n) with soft photon 
exponentiation: 

E^(m,n) = 1 + S{e, fi e ) + Hl{m,n) + Hf{m,n), (56) 

B^(m, n) = [1 + S(/3 e )]HI(m, n) + Hf (m, n). (57) 

The soft photon corrections »S(e, (3 e ) and S(f5 e ) have been introduced above. The hard photon parts 
consist of the residual hard part Hf(m,n) and another one, .fff 3( m > n ) > which depends on the 
treatment of the soft photon corner of the hard photon phase space integral: 

Hf(m,n) = giKiiRn) - Ki(I%)], i = l,3. (58) 
Here, g is a resonating kinematic factor: 

(1 - Rm)(l - i?*) , . 

9 = N -R* ' (59) 

The threshold behaviour and further details of the dynamics are contained in the JC-functions: 

1 



K\(R) 

K 2 (R) 
K 3 (R) 
I{z) = - 



-8 e R[A+(l + R)I(l-R)}, 



+-. 



R 
R 



p e [RI(l-R) 
A 



mi 



R 



A^J 



1 



dv- 



1 



In 



(60) 

(61) 
(62) 
(63) 
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J(a) = 2 F 1 (l,p e ,l + f3 e ,a)=(3 e 



1 yPe 



dv 
o 1 



av 



(64) 



The above definitions are valid if at least one of the R m , R n are non-zero. 

Although the limit of both masses M m , M n vanishing formally exists, it differs from the correct 
form of R^(0, 0) by a constant. This is due to the occurence of the double pole (1 — v)~ 2 from 
s'; 0, 0) under the integral. Such a double pole behaviour cannot be calculated as continuous 
limit of the two single poles arising in the massive case: 



i*r(0,0) 




~ S ] + -QliLe 
7T 



i) 



^2/ 1 



A' 



(65) 



where the upper, barred (lower, unbarred) case is with (without) soft photon exponentiation. Fur- 
ther, if the cut-off A is removed, an additional finite error would occur due to lacking phase space 
factors ( ^). We quote here the corresponding correction without E) and with soft photon 
exponentiation: 



R e T (o,o) 




1 + S(f3 e ) + 52 +S ] + —Q 2 e (L e 
IT 



1) 



1\. s 2/1 



p = 1 — 4mj/s. 

The logarithmic final state mass dependence regularises the massless photon propagator for emis- 
sion of very hard photons and the corresponding singularity of ln(l-A) in the limit A — > 1. For 
applications, we add to the soft photon correction S{(3 e ) the next to leading order logarithmic cor- 
rection 5% +S exactly as obtained in [|H>[Hl. For Hf(l,l) which is the dominating hard photon 
contribution at LEP energies, we also add terms aris ing from an explicit integration with the corre- 
sponding hard photon terms in 5^ as quoted in 113 113 . This ensures the necessary high accuracy 
of the order 0(0.1%) but is not described here. 

Due to the occurence of logarithms in the hard radiator function for <jfb (see 113 for details), 
the corresponding integrated expressions have a slightly more complicated structure than those for 
ot • For the photon exchange contribution, 

R e FB (0, 0) = i%(0, 0) + ^Q 2 e {~(L e - 1) In ^ 

+4d + (L e - l)[Li + L 2 + 2d] + 2 j^Li + ^L 2 + k + / 3 j . (66) 

The following abbreviations are used: 

d=lw^ir, (67) 
22-A' v ' 

Li=ln(l-A), (68) 
L 2 = In (l - I) , (69) 

h = Lia(A) - 2Li 2 (j) , (70) 
1\ ../IN . „ La 



*4 = Li 2 ( -J - Li 2 f—j - ln2L 2 + -L 2 , (71) 
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P 1 rlr 

Li 2 {z) = - / — hdl-xz). (72) 
Jo x 

The initial state correction to ctfb due to the 7Z interference is: 

*™a,o)=(*f ){ j[A/(i r i?)] } [ i + ^ e)+ ^] 

1 ) (^-l)[lnA + L iJ (A)] 
+2(L e - l)[d r - ^L R (A)} + 4 d r + 2r]^-^L 1 



+ 



(L ' l) (1 + 37?) + 4 d 



rL 2 



+(1 - 2i?)/ 4 + Z 3 + 2i? , D , 2 ^3 • (73) 



1 + 

1 + i? 2 

The resonance logarithm L R (A) as defined in ( ^5|) and the additional abbreviations depend on the 
complex Z boson mass parameter: 

D 3 = 2D 2 -D l + U + (L e -l-2 In 2)L R (A), (74) 

B2 - Lh {ttr) - Lh (it!) + ln21n( - r) - (i2 + ln 2)ln (' - frl) • (T5) 

z^^I)-^)-^-!^), ,ra, 



1-R 

f = — . 

1 + R 

The third correction to ctfb is due to Z exchange: 

A^ e 



(77) 



V+Si 

2 J 



2a o2 I / 



+— x (L e -l)[lnA + tL i? (A)] 



-t (L e - 1)Lb(A) + 2|rfd(L e + 1) + 2R .] + f. 2 tD 3 + 2|r| 2 ^J^ jLl + Ad\r\ 2 L 2 



+(L e - l)L 2 [8|^-^f + + 5] + (-1 + 2|1 - i?| 2 )/ 4 + h}, (78) 

The C-odd initial state corrections behave very similar to those for cry . They have also the loga- 
rithmic electron mass singularity L e . Both Z exchange corrections R%(1, 1) develop a radiative tail 
beyond the resonace. This is due to their structure as already discussed for the angular distribution; 
see ( ||): 

R A (1, 1) = f\{R) + fl(R)tL R (A)(L e - 1). (80) 

Because of their complexity, we quote here the functions R A (1,0) and R A (1, 1) for A = 1 taken at 
R = 1, i.e. at resonance (we leave out here the coupling constants in the definition of x ( HH)) : 



\ X \ 2 Re R e A (l, 1)\ R=1 ~ 4[1 + ?Qe(M + + Aa)], A = T, FB, (81) 



7T 
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P = T z /M z , (82) 

K = -\ + 2Li 2 (1) + (L e -l)[| + 21np], (83) 

h e 2 = -27r(L e -l), (84) 

fc§,= (L e -l)(2-3bip), (85) 

h FB = \ ~ + 4 In 2 2 - ^Li 2 (l) - i hip + (L e - l)(-~ In 2 + I - | hip). (86) 

From these 0(a) expressions it is easy to derive the corresponding ones for soft photon exponenti- 
ation. The most interesting feature is connected with hard photons - or better, with their absence 
at resonance. The initial state emission of a hard photon for s = M§ leads to a largely reduced 
effective energy s' and thus a non-resonant behaviour, i.e. a much reduced cross section. As a 
consequence of the resulting soft photon dominance, C-even and C-odd observables behave similar 
and consequently the leading order coefficients h\,h\ are equal for them at resonance. This has 
been observed numerically first in @3 and explained in —I. The 7Z interference corrections at the 
peak position are exclusively due to the imaginary parts of R e A {l,Q) since x becomes imaginary for 
R = l: 

Re[ X R e A (l, 0)] ~ -{5 A>T + -Q 2 ML e - 1) + pg\ + p 2 g A ]}, (87) 

P 7T 

<# = 2(L e -l)(lnp-i), ^ = _|(L e -l), (88) 
g° FB = 2{L e - 1) (In p - - + In 2) + [-2 + 2 In 2 - 4 In 2 2 + Li 2 (1)] , (89) 

9fb = ~J L e- (90) 

The Born contribution to &FB vanishes at the peak. Nevertheless, the leading corrections to o^and 
(Jfb are equal again. Integrated cross section and asy mmetry as functions of energy and the photon 



energy cut-off A have been studied in detail in O E3 . We will come back to that point in chapters 
5 and 6. 



4. INITIAL-FINAL STATE INTERFERENCE RADIATION 



The initial-final state interference corrections have an interesting property: only those which 
are diagonal in the arguments indicating the exchanged vector bosons are independent observables. 
The other interference contributions may be determined by the following simple relation: 

R A (c;m,n) = ^[R A (c;m,m) + ii^(c;n,n)*] (91) 

In order to get ( |9l]), it is essential to separate the reduced Born factor ( 1C, ^) from the QED 
contents. The validity of ( |9l| ) may be seen immediately from ( |3[ |l2[ |l3|). An explicit proof has 
been given in Hi! . 



4.1. THE RADIATOR FUNCTIONS FOR THE ANGULAR DISTRIBUTIONS 
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For interference corrections, the soft photon part has to include besides the soft photon emission 
terms Sa{c, e, A) also the contributions B A (c, A; m, n) originating from photon - Z boson and photon 
- photon box diagrams. Their sum is infrared finite (and independent of the infrared cut-off A) while 
the infinitesimal soft photon cut-off parameter e disappears only after integrating over the photon 
energy variable v: 

OL 

S l A (c, e; m, n) = -QeQf[SA{c, e, A) + B A (c, A; m, n)j, (92) 
S A (c, e, A) = 2D A (c)[2\n - In — + Li 2 (c+) - Li 2 (c_) - ^(ln 2 c+ - In 2 c_)]. (93) 



Here, £|(c) is defined in ( g|- ||) with A = T(FB) il A = FB(T), and c± in ( gT]). The box 
diagram contributions M are dependent on the type of the gauge bosons exchanged. We write 
them in (anti-) symmetrised form: 



Bt,fb(c, A; m, n) = b(c, A; m, n) ± b(—c, A; m, n). 
The two different box functions are: 



b(c, A; 0,0) 



-2c 2 In ^[4 In ^ 

c_ A 



2m] + lnc_[— 2c_|_ + c(lnc_ + 2m)] — 2nic + , 



(94) 



(95) 



b(c, A; n, n) 



1 ( 



= -4c+(l - R n ){\n — - (1 - i2 n )L n + 
1 - i?„ - 2c+)[/(l) - Z(c_) - L n lnc_]} 



c 

-2c 2 {[2 In — + AL n + m(c + c_)] In -± + 2Z(c+) - 2Z(c_)}, 
A c 

The following abbreviations are used: 

1(a) = Li 2 (l — aR~ l ), 
L n = ln(l-R- 1 ) = L Rn (l), 



and R n is defined in ( L R (A) in ( [45|). 

The hard radiator parts are independent of the gauge boson exchanged: 



a 



h t,fb( v , c ) = -QeQf[h % T FB (v,c) ±h l TFB (v,-c)}. 

7T 



(96) 

(97) 
(98) 

(99) 



While the box terms for A = FB and A = T could be expressed by one and the same function, 
this is not the case here: 



hir(v,c) 
2c+{ 

+2 



3- z(2 + z) 



-f + 4 + z(2 + lnz) 



l n £^_(l + z )2 ln £=±££+\ 
2-2-3(1 + *) 



+ ^ 

7 



(100) 



h F B(v,c) = 

2(l + c 2 )lnc_(f -1-z-z 2 ) 
+4c + [-4 + 4 + 23 - 3(1 - z) In*] + + § - 3 - z) 



+ | " 4 - 2z - z 2 ) + 2(z 2 - l)c + In c+c_ + 



2[(l-z + z 2 ) + c(l 



+ c 2 (l + 3 + z 2 )] In 7. 



(101) 
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It is easy to check that in the soft photon limit (v — ► 0) the functions h A (v, c) behave such that they 
compensate after integration over the photon momentum the dependence of the box contributions 
on lne. In Fig. 3, the interference contributions are shown as functions of the scattering angle for 
three different energies. Compared to initial state radiation, the interference corrections are small. 
They do not contain fermion mass singularities though kinematic singularities at cosO = ± 1 occur. 
In contrast to initial state radiation, these are integrable so that formally the full angular range in 
( p^, 24) may be used. Of course, the interference contributions depend on the hard photon cut-off 



A. For more severe cuts, they grow up. This is due to the fact that for the totally inclusive problem 
(A = v m ~ 1) there exists a fine-tuned cancellation of box and bremsstrahlung contributions (for 
more details see next section) which becomes more and more disbalanced if the phase space of the 
bremsstrahlung integral becomes more and more restricted. For very tight cut values, A <C 1, the 
cross section to order 0(a) starts to diverge and may become even negative. 

This can be seen from the leading soft photon contribution i?^ so '*(c; m, n). After integration 
over the photon energy, the dependence of S A (c, e, A) in ( [9^) on e drops out. It is replaced by the 
following terms arising from the soft photon corner of the phase space integral: 

fi%° f \cim,n) = 2D A (c)f3ilnA, (102) 

QJ 1 — Q 

P i = 2-Q e Q f \n-—. (103) 

7T 1 + C 

The complete interference contributions to the differential cross section may be found in the code 
MUCUTCOS In analogy to initial state radiation, the logarithmic dependence on the cut- 

off A in ( |102| ) has to be cured by an adequate soft photon resummation which should lead to a 
replacement of 

/ dv afis^s'^m^mv^pdne + Qiv-e)— 1 (104) 
Jo v 

by some function of the following type: 

dm<TT{,8,sf',m,n)PivP i - 1 . (105) 

o 



In order to get really a smooth, well-defined behaviour of ( |105 ), the exponent therein has to be 



larger than -1, i.e. fa should be positive definite. This is not the case in ( |105| ) as may be seen 
from ( |103j ), A possible way out is the exponentiation of interference radiation together with initial 
and final state radiation as is dictated by soft photon theorems (see, e.g., 0). A collection of the 
relevant terms leads to the following order O(a) expressions after integration. (One should have in 
mind that s' = s for the soft photon case): 

rlrr so f t 9rv s 1 — r s 

^7" = E <(s,s ] m,n)D A (c){l + —[QlH^-l) + 2Q e Q f ln—^ + Q 2 Mn^-l)}lnA}. 
dc ~ 7r mi 1 + c J mi 

A,m,n e j 

(106) 

The final state expressions have been taken from sect. 5. Again, the correction is not positive definite 
due to the expression in square brackets wich is a possible candidate to replace the exponent in 
(|i~05D. This is a consequence of too crude approximations. The soft photon contribution as a whole 



must be positive for any parameter combination since it is a complete module squared. In fact: 

da so ^= £ a A (s,s;m,n)D A (c)(l + - X 2 ), (107) 



X 



\kip k 2 pj \pip P2VJ 
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The four momenta used here are defined in ( |I|). After a little algebra, one gets 



X 



-2ml 



Qe L(-*ip)(-*2j0 (-fcip) 2 
+QeQf{t 



+ Q 2 f 



-2m 



f 



1 , 1 

(-k 2 p){-pip) ~ {-k 1 p)(~p 2 p) 



(-piP)(-P2P) (-PiP) 2 



7 



(-P2P)' 



1 , 1 

(-felP)(-PlP) (-fc 2 p)(-P2P) 



Am 2 e Js — Am?j. 



(109) 
(110) 



s c 

+(-) = -2pi(2)*2 = -2p2(i)fci = g ± 2 
er integration over the soft photon corner of the n dimensional photon momentum phase space 



153, one gets for the corrections which are relevant for soft photon exponentiation: 



da soft = ^A(s,s;m,n)D A (c)ll + -^-Q 

A,m,n 



2a ^ „ 
+ — QeQf 
IT 




(T , cr, + R 
-2- In — 

R~ ffr - R 



(s 2 -Am 2 e(f) s) 



1/2 



±2 
c 



s 2 db ca e af, 

2„2^1/2 



16m e mjS" 



(111) 

(112) 
(113) 
(114) 

whole is positive definite, 



In the ultra-relativistic limit of ( 111), one gets ( |106j ). While ( |111 
this property is lost in ( |106j ) as a result of the approximations. 

Since we do not intend in our applications to apply such stringent cuts as to make inevitable 
the adequate exponentiation of initial final interference contributions we would like to remain at 
this point of clarification. 

In a recent series of papers, there have been presented quite interesting results on soft photon 
exponentiation eiJ including the interference ea. Tracing back also to earlier work on that subject 
the problem of exponentiation of soft photon interference radiation has been solved there in 
the ultra-relativistic limit. From the above discussion we conclude that this is adequate for values 
of cosQ not too close to 1, i.e. assuming some realistic acceptance cut for the experimental set-up. 
Further, the authors make the implicit assumption that one may use the radiators for ot also for 
the differential cross section. This is not the case in general as shown in this article and in 03. 
Nevertheless, the radiators for ut can be used as an approximation. Further, we see no problem in 
combining the hard radiator parts presented here for the differential cross section with the result 
of E3 for the refined treatment of soft photon exponentiation. 



4.2. INTEGRATED CROSS SECTION AND FORWARD BACKWARD ASYMMETRY 



The integrated interference contributions are composed of box corrections Ba(\, n, n) and of 
bremsstrahlung contributions 6^(A, n, n): 

a 

-Q e Q f [B A {\; n, n) + b A {\; n, n) . (115) 

7T 
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For the total cross section one gets: 



B T (X;0,0) = 6 In 



2E 9 
X" 2' 



2_E A 

MA; 0, 0) = -6 In — — + 3(2 + A), 
A 

2E 

B T {\- 1, 1) = 6 In — - 9 + 3R[1 + (1 + R)L R (1)\ + 3L Z , 
A 

2 FA 

6 T (A; 1, 1) = -6 In -— + 6 + 3[A(1 - R) - R(l + R)L R (A)], 
A 



L Z = ln( S /M^), 



(116) 

(117) 

(118) 

(119) 
(120) 



where Lr(A) is defined in ( 
to <tt are quite small H 1 : 



and M' in ( ^3|). Without cut (A = 1), the resulting corrections 

(121) 
(122) 



a 9 
RU0,0) = -Q e Q f -, 

TT Z 

a 



Wr(l,l) = -QeQfiLz. 

TT 

The interference QED corrections due to photon exchange to the antisymmetric cross section part 
@fb are: 

2E 3 1 
Bj. B (A;0,0) = (1 + 8 In 2) In — - -(1 + 6 In 2 + In 2 2) + -i?r(2 - 5 In 2), (123) 

A 4 2 



IE A. 1 1 
Vf B (\;0,0) = -(1 + 81n2)ln— — + -(2 - 5A) + -(4 + 26A - 5A 2 )ln2 

A 4 4 



i 



i 



+-(-44 + 26A - 5A Z )L 2 + -(15 - 20A + 5A 2 )L X 
4^2(1) + 41n 2 2 - 4Li 2 (A) + 8Li 2 ( —) - -Lt 2 (A - 1). 



Their sum is a simple constant for A = 1 c3 : 



a 



R l FB (0, 0) = -Q e Q/(4.572 - »2.302). 



7T 



(124) 



(125) 



The Z exchange corrections are much more involved. For A = 1, they may be found in with 
A < 1: 

IE 3 

5^ B (A; 1, 1) = (1 + 8 In 2) In — - - + R - (9 - 4R - 4R 2 ) In 2 

A 2 

-2 In 2 2 + -(5 - 4R)L Z - - [4 - 9R + 3i? 2 + 2(-5 + 3i? - 6i? 2 ) In 2]L R (l) 



1 



1 



-(1 -3R + 6R 2 - 8iT) [Li 2 (l -) - Li 2 (l - — 



2R' 



R 



4iT 



Li 2 (l)-Li 2 (l-i) 



(126) 



IE A 1 

6 J Fi3 (A; 1, 1) = -(1 + 8 In 2) In — — + -(2 - 5A + 5i?A) 

A 4 

+^(4- 16i?A + 5i?A 2 - 10 J R 2 A + 26A-5A 2 )ln2 + 41n 2 2 + 2Li 2 (l) 
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+|(3 - R + 2RA - RA 2 - 2R 2 + 2R 2 A - 4A + A 2 )Li 

+(_ 5 + 3R _ 4i?A + ^ftA 2 + 5R 2 - b -R 2 A - -A- + H A - ^A 2 )L 2 
4 2 1 + i? 2 4 

+fl(-5 + — & — + fl)Lji(A) - 4Li 2 (A) + 8Li 2 £) 
i + -ft i 

+^(3 + 5R 2 )D 1 --(5 + 3R + 3R 2 + 5R 3 )D 2 , (127) 



with D\ t 2 defined in ( [75|- 76). Again, the Z boson parameters have to be understood in primed 
quantities. A generalization of the initial final interference QED corrections to the case of several 
massive gauge bosons is trivial due to relation ( |9l| ). At the Z peak, a fine-tuned cancellation of 
box and bremsstrahlung terms occurs. Therefore, the resulting Z exchange interference contribu- 
tions Rrpjl, 1) and R l FB (l, 1) become small there. This has been noticed from numerical results in 
IB @' 13 anc [ j s known also from earlier investigations of 3/ip physics For 1) it is evident 

from ( |122| ). A Taylor expansion around s = M| yields: 

= -Q e Q/[61n - - \ + A + 9(R - 1) ln(R - 1 + A) + 0(R - 1)], (128) 

7T A 

R i FB (l, 1) = -Q e Q/{[l + 81n(2 - A) - 4m(l - A)] In - - \ + A - 3Qn2)(R - 1) ln(R - 1) 

7T A 

+-[1 + 81n(2 - A) - 61n(l - A)](R - 1) ]n(R - 1 + A) + 4(R - 1) ln(l - A) + 0(i? - 1)}. (129) 
For A approaching 1 and using in the peak region In R = R— 1 + 0[(R — l) 2 ], one gets 

i4(l, 1; A = 1) = ^^{-3(12 - 1) + 0[(R - l) 2 ]}, (130) 

i^ B (l,l;A = 1) = -Q e Q f {-3(R- l)(ln2) ln(i? - 1) 
+^[7- 101n2-61n 2 2-6Li 2 (l)]( J R- 1) +0[(fl- l) 2 ]}. (131) 

The real parts of i?^(l,l) contribute to a a in a product together with |x| 2 ~ |i? — 1| -2 - As a 
cosequence of their proportionality to {R — 1) whose real part vanishes at the peak, the Z ex- 
change interference corrections become extremely small. For A /l, there are at the peak larger 
contributions due to ln(i? - 1 + A) ~ In A in ( [j]|§ |j~29l) . 



A similar discussion applies to the 7Z interference corrections. Due to ( |9l]), half the R\(l, 1) 
are combined there with the resonating function x wich becomes imaginary at R = 1. So, the 
imaginary parts of -R^(l, 1) are relevant. Again, the QED corrections are suppressed compared to 
the Z exchange Born cross section behaving like |1 — R\~ 2 . 

So, all the interference contributions are small since the photonic corrections R A (0,0) are also 
non-resonanting at the peak. In case of cuts to the photon energy, A < 1, the 0(a) interference 
bremsstrahlung has to be taken into account properly. For very small A, logarithms of the type 
In A may become even dominating and then one should try to apply some soft photon summation 
procedure. 
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5. FINAL STATE RADIATION 



5.1. THE RADIATOR FUNCTIONS FOR THE ANGULAR DISTRIBUTION 

The final state radiation contributions to the differential cross section have a simple angular 
dependence compared to the expressions discussed in the preceding sections. The soft photon parts 
of the final state radiators S^(c,e;m,n) may be obtained from S{e,j3 e ) as defined in ( |34|) by 
replacing m e by the final state mass mf : 

S f A (c,e;m,n)=D A (c)S(e,(3 f ), A = T, FB, (132) 
where Da{c) is introduced in ( [35, |33|). The hard photon radiators are: 



H£(v,c) = -Q){D T {c)[H f {v) -3v] + Av}, (133) 



H f FB (v, c) = -Q}D FB (c)h f (v), (134) 



7T 

O 
— ( 

7T 

where the radiators Hf(v), hf(v) have been derived earlier for the integrated cross sections or and 
a FB 0: 

HAv) = 1 + (1 ~ V) [L f - 1 + ln(l - v)], (135) 
v 

h f (v) = - v)(L f - 1) + ln(l -v) + ~v 2 L f ]. (136) 

The integration of the symmetric radiator function (A = T) over the photon energy gives a contri- 
bution to the differential cross section: 

R^(c; m, n) = i(5 - 3c 2 ). (137) 

The corresponding anti-symmetric radiator integral vanishes: 

R f FB (c;m,n) = 0. (138) 

As has been pointed out e.g. in EH , the soft photon corrections ( |132j ) to even and odd cross 
section parts are equal and will, consequently, cancel for the forward backward asymmetry. After 



integration over the full photon phase space, eqs. ( |137| - 138 ) , there is minor influence on Ap B due 
to the correction in the denominator. Nevertheless, for small cut off values A <C 1, soft photon 
exponentation is recommended in order to get reliable numerical results. A naive addition of initial 



and final state radiation corrections as has been performed in 113 changes this dependence and leads 
to an overestimation of final state radiation due to an inbalance of corrections. In order to maintain 
the smallness of isolated final state corrections after combination with the large contributions from 
the initial state, one has to proceed more carefully. A treatment based on the intuitive picture 
of subsequent radiation of photons from initial and final states is proposed in subsection 5.3. An 
isolated soft photon exponentiation for final state radiation is analogue to the procedure described 
in 113 and consists of the replacements 

5(v)S f A (c, e; m, n) -» PfV^^DA^S^f), (139) 
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H{(v,c) 



H f A (v,c)=H f A (v,c)-^D A (c). 



(140) 



Numerical results due to final state radiation are shown in fig. 4. 



5.2. INTEGRATED CROSS SECTION AND FORWARD BACKWARD ASYMMETRY TO 
ORDER 0(a) 



The integrated final state 0(a) QED corrections are simple: 



R A (m,n) = —Q 

7T 



a 



(L f - l)(41n A + 3 - 4A + A 2 ) - 2Li 2 (A) + 2Li 2 (l) + r{{A) 



4(A) = ~(l-A)(3-A)L 1 + i(6A-2 



A 



^ B (A) = -(A 



1). 



If the photon energy cut is removed, one gets: 



R f T (m,n) = "<2/^, 



IT 



(141) 

(142) 
(143) 

(144) 



, FB (m,n) = 0. (145) 

Figs. 5 and 6 illustrate the above discussion of the corrections to the angular distribution due to 
final state radiation. A comparison of ( |138| ) with ( |136j , [140D shows that the contribution of final 
state radiation to cjfb vanishes if neither a cut on the photon energy is applied § nor finite mass 
effects are taken into account e3. Then, the only (and minor) influence of final state radiation 
on Afb is due to the C-even correction in ( ||). If a tight cut is applied (A <C 1), the C-even and 
C-odd final state corrections ( |141[ - |143j ) approach each other because from soft photons there is no 
influence on the angular behaviour of the emitted fermions. Due to the logarithmic cut dependence, 
one again should exponentiate the soft photon contributions in that case but we leave out the details 
here. The ansatz is defined in ( |139 - |140j ). 



5.3. SOFT PHOTON EXPONENTATION 



As may be seen from ( |l 33| - 1136| ), the final state radiation contributions to the differential cross 
section behave as a function of the normalised photon energy v quite similar to the integrated 
quantities cr A ,A = T,FB. Based on their simple dependence on the scattering angle, we now 
derive a common treatment of soft photon exponentiation for initial and final state corrections. We 
start from the definition ( |3|): 

d a (e+f) ( /-A 

— = X] Re \ / dv a A (s',s';m,n) 

aC m,n=0,lA=T,FB l J ° 

{5(v) [1 + S(e, P e ))D A (c)} + 9{v - e)H A (v, c)} 
+<T A (s,s;m,n)D A (c) f du[6(u)S(e, /?/) + 6{u - e)H f A (u,c)}\ , (146) 
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H{(u,c) = D A (c)H f A (u,c), (147) 
where H FB is independent of the scattering angle and H T has an almost negligible dependence. 



Further, ( [146[) can be rewritten as follows: 



dc 



i?e / (if s'; m, n) 

m,n=0,l A=T,FB 



+ S(e, &) + fj(s', c, A)]£> A (c) + 0(w - e)H%(v, c)}, (148) 

/•A 

r f A (s,c,A) = du[6(u)S(e,P f )+6(u-e)H f A (u,c)}. (149) 



o 

The dependence of f A on c is very weak for A = T and absent for A = FB. Now the following 
ansatz seems to include a quite reasonable approximation of higher order soft photon corrections: 

^!1= £ £ Re f A dva A (s f , S f ;m,n)F e A (v,c)F{(v,c,s'), (150) 

dC m,n=0,l A=T,FB J ° 

F%(v, c) = 5(v)D A (c)[l + S(e, &)] + 0(i> - e)H%(v, c), (151) 
F;J(7;,c,s') = l + f A (s',c,A). (152) 
The exponentiation of soft photon contributions leads to the following replacements: 

F A (v,c) - F|(«,c) = D A {c)p e vP<- 1 [l + S(l3 e )] + H%(v ) c), (153) 



A' 



f{(v,c,s')^f{(v,c,s') = J du^iA-^l + ^^+fl^^c, #)}, (154) 

In <t' 

P' f = — Q2(in — -1). (155) 
Here again the residual hard radiator parts H are 



H A (v, c, 0) = H a A (v, c, 0) - ^D A (c). (156) 



The integrations over the variable u can be performed in ( 154): 



F^v,c,s > )=A' f3 'f[l + S(P' f )]+G A (v,c, S / ), (157) 



G A (v,c,s') = -/ f A'(A' -4) - ^Q}[2Li 2 (A')+g A (A',c)], (158) 
g T (A', c) = i(l - A')(3 - A') ln(l - A) - \a'(A' - 6) - (s - ^) ^. (159) 

WA,c) = V, (160) 

In the above formulae, we introduced the variable A' which is connected with the reduced s' in the 
same way as A with s: 

A' = 1 - , (161) 



A = (162) 
1 — v 



21 



if one wants to obtain from the above expressions for the angular distribution those for errand Apg , 
it is helpful to start with the observation that the Pi are independent of the scattering angle with 
exclusion of g^(A',c). Since that dependence is minor, i.e. less than \%Q) = 0.1%Qj , one can 
neglect it. Then the final state radiation factor is completely independent of the scattering angle 
and the integration over c has to be done for F^(v,c) only: 

4 6+/) = E Re [ A dva° A (s',s f ;m,n)FX(v)F^v,s'), (163) 

m,n=0,l 

n(v) = lJ Q 1 dcF e A (v,c). (164) 

This integration yields the normal convolution representation for the total cross section a a, but 
now modified by the additional factor F A . So, by a direct Feynman diagram calculation we have 
shown that ( |163| ) is an approximation of sufficient accuracy though some well-defined terms of the 
order 0(a) have been neglected. A similar formula in the case of gt has been obtained by other 
methods in @1. The dependence of or and Afb on the cut-off parameter A is shown in figs. 5 and 6 
for different handlings of initial and final state corrections. For smaller A the cross section becomes 
smaller since the positive hard photon part is restricted then. In the tail region (figs. 5c, 6c), one 
observes a steeply falling cross section at values of A which cut away the radiative tail. This happens 
if the effective energy s' cannot reach the resonance energy, i.e. A < 1-M^ /s. Even for infinitesimal 
A the net negative bremsstrahlung corrections remain finite due to soft photon exponentiation as 
discussed in the text. In the figures, we applied the simplified but sufficiently accurate procedure 
explained in chapter 4.2. Since the corresponding soft photon corrections in units of the Born cross 
sections <jt and (Jfb are equal, the net final state correction to Afb vanishes for vanishing A. For 
sufficiently large cut-off values, the final state contributions also remain nearly negligible. In the 
intermediate region, they must be taken into account at least for precision measurements. 



6. DISCUSSION 



In the foregoing chapters, basic formulae for and characteristic features of the QED corrections 
to three different observables in fermion pair production from e + e~ annihilation have been presented 
separately. NumericaLresults have been produced with the codes MUCUT and MUCUTCOS of 
the package ZBIZON©. 

In fig. 7, the net QED corrections to the differential cross section are shown in combination 
with the weak loop corrections. The latter have been determined with the code DIZET which 
is also part of ZBIZON. For M z =91.1, m t =100.,M H =W0. (all masses in GeV), and a s =0.12 
we get in accordance with sin 2 ©iy= 0. 2314, r^=2. 477 GeV. The corrected muon decay 

constant, running QED coupling, and weak neutral vector and axial vector couplings to be used in 
(0,0, g§fH]) are ©: 

aF A {M 2 z ) = a(1.063 - i0.018), (165) 
G^p z {Ml) = 6^(1.000 - i0.005), (166) 

«e = -~, (167) 
u e (Aff) = --[1 - 4sin 2 0^(1.011 + iO.013)], (168) 
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v ee (M 2 z ) = ~[1 -8sin 2 0^(1.011 + i0.013) + 16sin 4 0^(1.023 + i0.027)]. (169) 

If, for the sake of comparison, the initial- final interference terms are excluded, the fully corrected 
total cross se ctio n ot as obtained here agrees within ±0.2% with that obtained from the code 
ZSHAPE laila which is exact to order 0(a 2 ). For Afb and dcr/dcos© , we do not know of a 
program which would allow for a similar comparison^ 

While the weak loop corrections remain small 113 US, the QED corrections amount typically 
to several percents or more, essentially due to initial state radiation. Even in the tail region, the 
angular dependence of the corrected differential cross section follows closely the Born cross section 
behaviour with exclusion of scattering angles near |cos9| = 1. This is reflected also in figs. 8 
and 9 where ot and Aps are shown as functions of an acceptance cut, |9| < Q ma x- In general, 
the asymmetry is not a monotonic, rising with Q max function though there are regions with an 
almost linear behaviour. For loose cut values, the tail effect is pronounced at s > M§. The total 
asymmetry is small at resonance with corrections nearly of the order of the Born contribution. The 
relatively large asymmetry value for a tight cut, E~ <lGeV, may be understood from figs. 2b where 
a systematic distortion of the differential cross section is present in that case. 

In some recent determinations of ot 8 and Afb ^3 at LEP, the QED corrections '11' HBUll 
have been applied to data corrected to correspond to a detector with a full angular acceptance (Air 
geometry). From figs. 8,9 one can see that this may be justified although the present formulae 
allow a refined analytic approach to the data. This seems to be recommended at least for the 
interpretation of precision measurements. 

In view of the relative smallness of interference and final state corrections, it is interesting to 
know where they can be neglected completely. Of course, this depends on the energy region and 
photon energy cut-off chosen. In figs. 10,11 corresponding regions of relevance are shown. 

To summarise, the analytic approach to QED corrections has been proven quite powerful both 
for consistency checks of Monte Carlo programs H and for the interpretation of experimental data. 
This article contains a first, systematic presentation of our analytic formulae for the complete set 
of QED corrections to the differential cross section dcr/dcos©, total cross section errand integrated 
forward backward asymmetry Afb ■ 
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FIGURE CAPTIONS 



Figure caption Fig. 1. The C-even and C-odd contributions to du/dcosO due to Z exchange 
(solid line), 7 exchange (dashed line), 7Z interference (dotted line) in nbarn as functions of the 
scattering angle at y/s =30 GeV (a) 91.1GeV (b), 200 GeV (c). Parameters: M z = 91.1 GeV, T z 
=2.5GeV,sin 2 Q w =.23, A=l. 

Fig. 2. The differential cross section with Born and initial state radiation contributions in nbarn 
as functions of the scattering angle and of the photon energy cut-off A. Other parameters: see fig.l. 

Fig.3. QED contribution to da/dcosO due to initial-final interference as function of the scat- 
tering angle and A. Parameters: see figs. 1,2. 

Fig. 4. Born plus final state correction contributions to da/dcosO in nbarn as function of the 
scattering angle and A. Parameters: see figs. 1,2. 

Fig. 5. The total cross section ax as function of the photon energy cut-off A. Parameters: see 
figs.1,2. 

Fig. 6. The integrated forward-backward asymmetry Afb as function of the photon energy 
cut-off A. Parameters: see figs.1,2. 

Fig. 7. The differential cross section with complete electroweak corrections as function of the 
scattering angle and A. Parameters M^=91.1, m t =100.,M#=100. (all masses in GeV), a s =0.12. 

Fig. 8. The total cross section <jt as function of an acceptance cut on the scattering angle, 
I cos ©I < cosO maa; . Parameters see fig. 7. 

Fig. 9. The integrated forward-backward asymmetry Afb as function of an acceptance cut on 
the scattering angle. Parameters see fig. 7. 

Fig. 10. Cross section contributions due to initial-final state interference (a) and final state 
radiation (b) in percent as functions of yfs and A. Parameters: see figs.1,2. 

Fig. 11. Asymmetry contributions due to initial-final state interference (a) and final state 
radiation (b). Parameters: see fig. 10. 
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